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I  EXECUTIVE  SUMMARY 


The  disruptive  effects  of  naturally  occurring  propagation  disturbances 
on  conmunlcation  and  navigation  systems  are  well  known.  More  recently, 
however,  antenna  aperture  synthesis  techniques  are  being  considered 
seriously  for  space-based  surveillance  and  ground-based  space  surveillance. 
The  applications  include  both  passive  and  adaptive  schemes.  We  shall  refer 
to  such  systems  as  Discrete  Focused  Arrays  (DFA). 

To  support  the  planning  and  evaluation  of  these  new  systems,  the 
Air  Force  Geophysics  Laboratory  (AFGL)  has  tasked  its  scintillation 
research  program  to  address  the  impact  of  naturally  occurring  propagation 
disturbances  on  such  systems.  This  report  draws  on  many  facets  of  ongoing 
scintillation  research  to  identify  the  critical  parameters  for  the  problem 
and  to  review  the  recent  data  that  have  been  obtained  to  quantitatively 
define  these  parameters. 

The  propagation  theory  is  now  sufficiently  well  understood  that  we 
can  characterize  the  signal  structure  fairly  accurately  in  terms  of  a 
statistical  description  of  the  in  situ  irregularity  structure.  The  in 
situ  spectral  density  function  is  a  sufficient  statistic.  In  fact,  for 
most  engineering  and  many  scientific  applications,  the  ionospheric  distur¬ 
bance  can  be  replaced  by  an  equivalent  phase  perturbation  without  loss  of 
generality. 

The  DFA  problem  can  be  separated  into  two  parts.  First  the  propaga¬ 
tion  theory  is  applied  to  determine  the  disruptive  effects  in  the  specific 
systems  of  interest.  We  have  used  a  model  (prototype)  that  encompasses 
the  essential  characteristics  of  any  specific  system.  The  propagation 
effects  depend  functionally  on  a  mathematical  description  of  the  irregu¬ 
larity  spectral  density  function,  the  most  important  parameter  being  the 
perturbation  strength. 

Second,  the  conditions  under  which  severe  performance  degradations 
actually  occur  must  be  determined.  This  is  the  most  difficult  task  because 
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it  demands  a  detailed  knowledge  of  the  irregularity  structure  and  its 
spatial  and  temporal  morphology.  In  the  report  we  first  describe  the 
DFA  model  and  review  the  relevant  parts  of  the  propagation  theory.  The 
degradation  of  DFA  performance  involves  two  aspects  of  the  signal  struc¬ 
ture.  The  temporal  variations  of  the  signal  caused  by  the  scintillation 
mainly  affect  the  sensitivity  of  the  DFA.  Fading,  for  example,  reduces 
target  detectability.  Similarly,  phase  scintillation  limits  the  useful 
coherent  integration  time. 

The  synthesized  antenna  beam  is  affected  directly  by  the  spatial 
coherence  of  the  wave  field.  The  critical  parameter  is  the  spatial 
coherence  scale  compared  to  the  effective  size  of  the  synthesized  aper¬ 
ture;  however,  the  spatial  coherence  of  the  wave  field  depends  strongly 
upon  the  anisotropy  of  the  irregularities.  Our  calculations  show  that 
an  L-band  array  of  more  than  a  few  hundred  meters  extent  can  suffer  sig¬ 
nificant  performance  degradation  during  the  most  intense  disturbances. 

On  the  other  hand,  the  regions  in  which  such  high- 1  a t i tude  distur¬ 
bances  occur  are  highly  localized  in  latitude.  Static  models  tend  to 
exaggerate  the  likelihood  of  severe  disturbances  at  a  single  location, 
particularly  when  geometrical  effects  are  included.  Truly  predictive, 
dynamic  models  will  ultimately  be  needed.  Because  the  DNA  HILAT  satellite, 
scheduled  for  launch  in  summer  of  1983,  is  dedicated  to  this  problem,  we 
have  emphasized  the  scintillation  phenomena  that  are  relevant  to  global 
modeling  as  well  as  the  aspects  that  directly  relate  to  DFA  performance. 
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II  INTRODUCTION 


Distributed  antenna  systems  in  which  fixed  (nonsteerable )  elements 
independently  measure  phase-coherent  signals  and  transmit  them  to  a 
central  processor  are  being  considered  for  both  passive  and  active  sur¬ 
veillance  systems.  The  resolution  of  such  distributed  focused  arrays 
(DFA)  is  limited  by  the  beam  width  of  the  individual  elements,  because 
all  elements  must  illuminate  the  focus  point. 

To  obtain  a  simple  lower  bound  on  the  achievable  resolution,  let  L 

e 

be  the  effective  aperture  size  of  an  array  element.  The  maximum  element 
separation  that  can  contribute  to  a  focus  at  R  can  be  approximated  as 


e 


Similarly,  the  corresponding  spatial  resolution  of  the  DFA  is 

AR  =  =  L  (2) 

L  e 

max 

An  example  of  such  a  system  is  SEASAT  [ Jordan ,  1980],  which  uses  aper¬ 
ture  synthesis  with  a  satellite -borne  radar.  For  the  two-way  path,  AR  is 
half  the  effective  aperture.  SEASAT  can  achieve  a  AR  of  better  than  10  m 
with  digital  processing. 

For  distributed  focused  arrays  whose  focus  is  within  or  beyond  the 
ionosphere,  however,  the  harmful  effects  of  propagation  disturbances 
must  be  considered.  The  general  nature  of  propagation  effects  on  DFA 
performance  is  easily  understood.  Over  a  single  path,  the  propagation 
effects  multiply  the  complex  signal  and  cause  the  familiar  temporal, 
amplitude,  and  phase  scintillation  structure.  However,  because  the 
time-of-arrival  differences  to  the  various  DFA  elements  are  always  much 
smaller  than  the  scale  for  temporal  changes  In  the  scintillation  structure, 
only  the  spatial  structure  of  the  scintillation  affects  the  focusing  oper¬ 
ation. 
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The  array  is  focused  by  forming  a  weighted  summation  of  the  complex 
signals  at  each  element.  It  follows  that  any  scintillation-induced 
defocusing  can  be  described  by  a  convolution  in  the  spatial  wave-number 
(angle)  domain.  As  noted  above,  the  time-dependent  signal  variations 
do  not  affect  beam  formation,  but  they  do  limit  target  detectability. 

For  example,  fading  reduces  the  average  signal-to-noise  ratio.  Similarly, 
the  temporal  coherence  of  the  signal  limits  the  time  interval  for  effec¬ 
tive  coherent  integration  [ Rino  et  al .  ,  1977]. 

To  calculate  DFA  performance  degradation,  we  have  developed  a 
general  model  that  can  accommodate  any  DFA  configuration.  Propagation 
effects  are  described  on  each  path  by  using  the  Huygens-Fresnel  integral. 
An  equivalent  phase-screen  relates  the  initial  disturbance  to  the  in  situ 
irregularity  structure.  The  model  is  used  to  calculate  the  average  dis¬ 
tortion  of  the  DFA  beam.  For  this  application,  the  phase  structure 
function  is  a  sufficient  character ization  of  the  average  propagation 
disturbance.  The  critical  parameters  are  the  spatial  correlation  scale 
and  the  anisotropy  of  phase  perturbation. 

To  verify  the  model,  we  have  analyzed  a  large  data  base  of  spaced- 
receiver  measurements.  Spaced-receiver  measurements  are  important 
because  from  single-point  measurements,  neither  the  anisotropy  nor  the 
spatial  scale  of  the  wave  field  can  be  determined  unambiguously.  In 
fact,  there  is  a  complicated  interplay  between  the  anisotropy  of  the 
scintillation  diffraction  and  its  apparent  drift.  Spaced-receiver 
measurements  allow  us  to  determine  the  anisotropy  and  true  pattern 
drifts.  Temporal  correlation  scales  can  then  be  converted  to  spatial 
scales  and  the  effects  on  DFA  performance  accurately  predicted. 

In  addition  to  the  largo  body  of  spaced-receiver  data  that  have 
been  analyzed,  during  November  1980  we  performed  a  series  of  measurements 
using  two  independent  GPS  receiver  systems  with  antennas  separated  by 
100  m  along  an  east-west  baseline.  The  GPS  satellites  are  in  12-hr 
orbits.  Several  GPS  satellites  were  visible  at  high  elevation  angles 
for  many  hours  at  Poker  Flat,  Alaska,  where  the  measurements  were  made. 
Thus,  we  had  the  opportunity  to  search  for  scintillation  source  regions 
and  verify  the  DFA  model  for  a  single  pair  of  elements. 
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Ill  A  DISCRETE  FOCUSED-ARRAY  MODEL  FOR 
EVALUATING  PROPAGATION  EFFECTS 


The  model  described  in  this  section  was  developed  for  space-based 
radar  applications  [ Rino ,  1982a].  We  modeled  an  arbitrary  configuration 
of  noninteracting,  omnidirectional  elements.  For  purposes  of  analysis, 
it  is  convenient  to  think  of  the  elements  as  radiators.  The  nth  element 
emits  the  spherical  wave 


exp{-ikr  } 
n 


r 

n 


(3) 


By  reciprocity,  the  problem  is  invariant  to  a  source-receiver  inter¬ 
change  . 

If  the  nth  element  is  located  by  the  vector  A £  ,  then 

n 


r 

n 


aR(0,  *)  k 


(4) 


where 


(5) 


R  is  the  distance  from  the  center  of  the  reference  coordinate  system 
— *  —♦ 

to  the  tip  of  r  ,  and  aD  is  a  unit  vector  along  R.  The  angles  (0,  0) 
n  K 

— • 

locate  the  direction  of  R. 


We  shall  assume  that  6£  I  <  <  1  for  the  most  remote  DFA  elements, 

1  n 1 

Then,  the  approximation 


r 

n 


R  +  h 


aR(e,  +) 


(6) 


can  be  used.  The  array  is  "focused"  by  multiplying  each  element  by  the 

complex  weighting  factor  A  exp{i*  }.  The  DFA  field  distribution,  E„, 

n  n  K 
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is  computed  by  summing  the  contributions  from  each  element.  Thus, 


er(0,  +)  = 

^  (  r  I 1 ^  _  I 

*  Z  An  exp  |  i['*n  +  hk—fT~  -  kA?n  '  V6>  *>]  <7> 

The  quadratic  term  depends  upon  (0,  through  R,  but  only  weakly,  and 

can  be  absorbed  in  $  without  loss  of  generality.  We  can  also  ignore 

n 

the  spherical  reference  wave  and  consider  only  the  DFA  voltage  pattern 


VR(0, 


<*>)  = 


exp 


!-« I"** 

l  L  n 


kA£ 


V0 


■  *']! 


(8) 


The  array  is  focused  by  chosing  A#  to  cancel  the  linear  A £  term 

’  n  n 

for  some  fixed  direction  (0  ,  <t>  ).  In  that  case,  we  can  write  Eq.  (8) 

P  P 

in  a  simpler  form  as 

V(aK)  =  V*  A  exp  UA 1  *  AK}  (9) 

R  t—J  n  n 

n 

where 


AK  =  -k  aR(0p,  <0p)  -  aR(0 ,  <t>)  .  (10) 

Thus,  the  voltage  pattern  of  a  DFA  can  always  be  written  in  the  form  of 
a  discrete  Fourier  transform. 

To  accommodate  propagation  effects,  we  consider  a  phase  screen 
located  at  a  distance  H  from  the  array  elements  as  shown  in  Figure  1. 

The  propagation  disturbance  on  each  path  is  completely  characterized 
by  the  complex  time-varying  function  h  (t).  Under  even  the  most  dis¬ 
turbed  conditions,  however,  the  temporal  variations  of  h^(t)  are  very 
slow  (greater  than  tens  of  milliseconds)  compared  to  propagation  delays 
within  the  array  (less  than  tens  of  nanoseconds).  Thus,  insofar  as 
array  focusing  is  concerned,  temporal  changes  are  completely  negligible, 
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FIGURE  1  DNA  COORDINATE  SYSTEM  FOR  ANALYSIS  OF  DISTRIBUTED  SYSTEM 


and  we  can  write 

V  (AK)  =  V''  Ah  exp{iA«  *  AK}  .  (11) 

K  n  n  n 

n 

Even  so,  the  temporal  structure  of  the  scintillation  is  ultimately 
of  concern,  because  it  limits  the  sensitivity  (as  opposed  to  spatial 
resolution)  of  the  system.  The  fading  and  temporal  coherence  properties 
of  scintillation  signals  are  well  known,  however,  and  we  shall  concentrate 
mainly  on  the  limitations  to  resolution.  The  fact  that  this  partitioning 
can  be  performed  greatly  simplifies  the  problem. 

To  characterize  h^,  we  postulate  a  pure  phase  perturbation  at  H  and 
then  use  the  Huygens-Fresnel  integral  to  calculate  the  evolution  of  the 
phase  and  amplitude  structure.  Using  the  phase-screen  model  is  not 
unduly  restrictive.  Indeed,  the  average  scintillation  structure  of  the 
Wideband  satellite  data  can  be  accurately  reproduced  by  calculations 
based  on  the  phase-screen  model  [ Rlno ,  1982b J. 

The  integral  becomes  unwieldy,  however,  if  we  fully  accommodate  the 
angle  dependence  of  the  ray  path  to  each  element.  On  the  other  hand. 
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under  essentially  the  same  conditions  tliat  Eq .  (6)  applies,  we  can 
ignore  the  angle  dependence  of  the  individual  ray  paths.  If  the  array 
is  focused  on  axis, 


hn  K)  exp  J-iK2  ^  J 

,  .  -  .  dK 

x  exp  t-iK  *  p  C}  - 

n  (2* 

where 

C  =  1  -  H/R 

r  =  (1  -  H/R ) H 
e 

p  =  A?  (1  -  H/R) 
n  n 


(12) 


(13) 

(14) 

(15) 


and  u(K)  is  the  Fourier  transform  of  the  phase  perturbation  in  the 
plane  at  H.  The  general  expression  is  given  in  Rino  Ll982a], 

To  accommodate  the  temporal  structure  of  the  scintillation,  it  is 
only  necessary  to  replace  in  Eq.  (5)  by  A£^  -  v^t,  where  v^  is  the 

H-plane  component  of  the  apparent  irregularity  drift.  The  apparent 
drift  accommodates  antenna  and  target  motion  as  well  as  the  true  irregu¬ 
larity  drift.  We  assume  that  any  change  in  the  structure  itself  is 
negligible  over  time  intervals  of  interest.  With  these  assumptions,  the 
model  includes  all  potentially  deleterious  effects  on  DFA  operation. 

To  make  the  model  less  abstract,  we  shall  compute  the  average  DFA 

power  variation  or  gain  pattern.  For  typical  DFA  applications,  we 

attempt  to  direct  the  main  beam  and,  in  some  applications,  steer  the  side 

lobes  to  cancel  out  an  interfering  signal.  In  the  absence  of  propagation 

effects,  the  power  pattern  is  given  by  |V^(AK)|  . 

uniformly-spaced  linear  DFA  with  A  =1, 

n  ’ 

|  V  f  AK>  |  2  . 

s  in^  (AjBAK/2  ) 
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As  an  example,  for  a 


(16) 


The  beam  width  to  first  null  is  — ^  ,  or  in  radians  X/(2L).  The  main 

-_2  NAf 

side  lobes  fall  of  as  AK 

Because  the  propagation  disturbances  are  random,  we  consider  the 

^  2 

expectation  value  of  |V^(AK)|  •  By  direct  computation  from  Eqs.  (11) 
and  (12), 


< |  V  (AK)  |  >  =  h*>  A  A  exp|iAp  Ak| 

1  R  1  is  n  m  n  m  I  nm  1 


n  m 


where 


// 


<h  h*>  =  //  $(K)  exp] 
n  m  1  *  1 


|  J  dK 

{  -  V  K  Ap  C  ; 


nm  f(2x)2 


(17) 


(18) 


and  #(K)  is  the  spectral  density  function  of  the  initial  perturbation. 
In  terms  of  the  phase  structure  function, 


<h  h*>  =  exp 
n  m 


(19) 


where 


Dfi^(Ap)  =<[6$(p)  -  6<*>(p  ')]2^  .  (20) 

It  follows  from  Eq.  (17)  that  the  average  distortion  of  the  DFA  beam 
is  independent  of  the  Fresnel  radius.  The  only  dependence  on  the  propa¬ 
gation  distance  is  contained  in  C  [see  Eq.  (13)].  If  H  =  R,  then  C  =  0 
and  there  is  no  distortion  as  we  should  expect.  If  R  >  >  H,  that  is,  the 
focus  approaches  infinity,  then  C  =  1  and  the  propagation  effects  are 
maximized . 

To  carry  the  problem  further,  we  identify  the  Ap  contour  such  that 

D.  .(Ap)  =  1.  As  long  as  the  DFA  elements  lie  well  within  this  contour, 
o<p 

the  DFA  beam  distortion  will  be  small.  To  evaluate  this  contour  explicitly, 
we  use  the  approximation: 

D^(Ap)  2  G  f(Ap)mintp'1,  23  (21) 
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where  p  is  the  phase  power-law  spectral  index,  and  G  is  a  geometrical 
factor  that  depends  on  the  propagation  angle  relative  to  the  magnetic 
field. 


C?  «  r2X2£  C.F(p)  ,  (22) 

6<p  e  p  S 

and 

2  2 

CAp  -  B  Ap  Ap  +  AAp 

fZ(Ap)  =  - ^ '  (23) 

AC  -  B  /4 

In  Eq.  (22)  r^  is  the  classical  electron  radius,  X  is  the  wavelength, 

£  is  the  length  of  the  propagation  path  within  the  disturbed  region,  C 
P  15 

is  the  equivalent  isotropic  turbulence  strength.  The  complete  form  of 

Eq.  (22)  and  a  definition  of  f(p)  and  the  coefficients  A,  B,  and  C, 

which  depend  on  the  irregularity  anisotropy  and  propagation  geometry, 

can  be  found  in  Rino  [1982b]. 

The  important  points  to  note  here  are  that  the  average  DFA  beam  dis- 

2 

tortion  depends  only  on  the  perturbation  strength,  through  and  the 

— ♦ 

irregularity  anisotropy  through  f(Ap).  The  perturbation  strength  can  be 
determined  from  phase  scintillation  measurements.  The  anisotropy  can  be 
determined  unambiguously  only  by  using  spaced-receiving  measurements 
[Rino  and  Livingston,  1982], 


IV  POKER  FLAT  GPS  MEASUREMENTS 

A.  General  Morphology 

During  November  1980,  a  series  of  measurements  were  made  at  Poker 
Flat,  Alaska,  using  the  GPS  Navstar  satellites  as  signal  sources.  The 
satellites  are  in  12-hour  circular  orbits  and  they  transmit  broadband 
(PRN)  signals  at  two  L-band  frequencies,  1575.42  MHz  (LI)  and  1227.6  MHz 
(L2).  With  appropriate  receiving  equipment,  it  is  possible  to  measure 
amplitude  scintillation,  the  LI  to  L2  differential  phase  scintillation, 
and  the  total  electron  content  (TEC),  as  described  in  detail  in  Klobuchar 
[1981]. 

The  350-km  ionospheric  penetration  points  for  the  satellites  visible 
through  the  local  evening  sector  are  plotted  in  Figure  2.  The  satellites 
progress  from  the  north,  and  they  are  visible  for  more  than  two  hours. 

The  points  of  closest  approach  lie  to  the  west  of  the  Poker  Flat  observing 
station.  The  remaining  passes,  shown  in  Figure  3,  provide  more  extended 
coverage  but  they  occur  well  into  the  morning  sector. 

As  discussed  in  Section  II,  two  factors  contribute  to  DFA  beam  dis¬ 
tortion;  namely,  perturbation  strength  and  irregularity  anisotropy.  The 
GPS  observations  were  made  to  map  the  distribution  of  auroral-zone  irreg¬ 
ularities  over  large  latitude  and  local-time  sectors.  Most  other 
auroral-zone  scintillation  measurements  have  been  made  using  either  geo¬ 
synchronous  satellites,  which  are  observed  at  low  elevation  angles,  or 
low-orbiting  satellites  such  as  Wideband  and  TRANSIT.  The  former  give 
good  local  time  coverage,  but  poor  localization  of  structure.  Low- 
orbiting  satellites  give  excellent  latitudinal  coverage,  but  only  a 
coarse  local  time  sampling.  With  several  satellites  in  high-altitude 
orbits,  the  GPS  constellation  is  ideally  suited  for  nearly  simultaneous 
local-time  and  latitudinal  coverage. 
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LONGITUDE  —  degrees 

FIGURE  2  EVENING-SECTOR  GPS  PASSES 


The  most  recent  work  on  high-latitude  scintillation  morphology  has 
been  heavily  influenced  by  data  from  the  low-orbiting  Wideband  satellite 
r Fremouw  and  Lansinger,  1981;  Rasu  e  t  a  1 . ,  1981],  In  more  than  40  per¬ 
cent  of  the  nighttime  passes  recorded  at  Poker  Flat,  for  example,  a 
localized  scintillation  enhancement  was  observed  at  the  point  where  the 
propagation  path  lay  within  the  local  L-shell.  This  has  been  attributed 
to  a  geometrical  enhancement  from  sheet! ike  irregularities,  although  a 
localized  source  is  also  involved  |\soe  Rino  and  Vickrey,  1982,  for  a 
detailed  discussion].  In  any  case,  static  (average)  morphological  model 
necessarily  are  heavily  influenced  hy  this  geometrical  factor. 


LONGITUDE  —  degrees 

FIGURE  3  MORNING-SECTOR  GPS  PASSES 


Because  the  GPS  satellite  moves  very  slowly  through  regions  of  higher 
geometrical  enhancement,  we  expected  a  preponderance  of  enhanced  scintil¬ 
lation  events  in  these  regions.  A  very  different  pattern  was  actually 
observed.  We  maintained  a  reasonably  consistent  schedule  of  operations 
for  nine  consecutive  evenings.  After  that,  the  operation  sequence  was 
changed  to  accommodate  dual-receiver  measurements,  which  are  described 
in  Section  III. 

The  measured  LI  intensity  scintillation  for  the  first  nine  days  is 
summarized  in  Figure  4.  Local  midnight  at  Poker  Flat  occurs  at  10:00  UT. 
Superimposed  on  the  data  for  Day  315  is  the  geometrical  enhancement  factor 
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for  10:10:1  sheetllke  irregularities  [Livingston  et  al.,  1982].  Tlie  di¬ 
urnal  progression  of  the  satellite  passes  is  small  enough  that  essentially 
the  same  pattern  can  be  applied  to  any  evening. 

Hie  most  prominent  sustained  enhancement  that  can  be  associated  with 
geometry  occurs  on  the  GPS  4  path  just  before  800  UT.  A  small,  but  con¬ 
sistent  enhancement  occurs  on  the  GPS  5  path  around  0930  UT.  The  dominant 
scintillation  enhancement  occurs  near  1100  UT,  however,  where  the  geomet¬ 
rical  enhancement  factor  is  small.  The  persistence  of  this  feature  night 
after  night  clearly  suggests  that  it  is  caused  by  a  localized  irregularity 
source  region. 

It  is  convenient  to  plot  the  data  on  an  invariant  magnetic  latitude- 
local  magnetic  time  plot.  Figure  5  summarizes  the  first  eight  days  of 


FIGURE  5  LOCATION  OF  SCINTILLATION  ENHANCEMENTS 


operation  in  this  manner.  The  activity  increases  as  we  move  equatorward 
and  toward  midnight;  moreover,  the  region  of  the  most  prominent  scintilla¬ 
tion  enhancement  occurs  near  the  expected  location  of  the  diffuse  auroral 
boundary. 
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The  Chatanika  radar  has  provided  unambiguous  identification  of  local¬ 
ized  scintillation  enhancements  that  are  most  prominent  when  viewed  along 
the  local  L-shell.  This  has  been  interpreted  as  a  geometrical  effect 
from  sheetlike  irregularities.  On  the  other  hand,  the  most  recent 
Chatanika  studies  indicate  that  these  structures  are  highly  dynamic, 
possibly  even  transient  factures.  In  addition,  a  much  more  stable 
"boundary  blob"  or  wall  of  F-region  ionization  is  located  at  the  equator- 
ward  boundary  of  the  diffuse  auroral  region. 

The  location  of  the  persistent  enhancement  seen  on  GPS  5  is  consis¬ 
tent  with  the  expected  location  of  the  boundary  blob  as  predicted  by 
Tsunoda  and  Vickrey  fl982].  These  results  are  preliminary,  but  it  is  clear 
that  accurate  prediction  of  scintillation  enhancements  ultimately  will 
require  an  understanding  of  the  evolution  and  dynamics  of  these  source 
regions . 

The  activity  is  summarized  in  Table  1,  which  lists  the  fractional 
number  of  data  points  with  at  LI  greater  than  0.1,  the  maximum  value 
of  at  LI  for  the  pass,  and  the  College,  Alaska  3-hour  K  index  of  mag¬ 
netic  activity.  This  comparatively  small  data  set  did  not  reveal  a 
systematic  dependence  on  the  magnetic  activity  at  the  receiving  station 
as  was  found  in  the  Wideband  satellite  data  [Rino  and  Matthews,  1980], 
probably  because  the  22,000-km  height  of  the  satellite  generally  places 
the  ionospheric  penetration  point  well  removed  from  the  vicinity  of  the 
station. 

In  summary  then,  over  the  eight-day  period  where  a  consistent  GPS 
observing  pattern  was  maintained,  very  similar  occurrence  patterns  of 
enhanced  scintillation  were  observed.  A  small  enhancement  was  observed 
on  CPS  4  before  0800  UT,  which  is  consistent  with  the  expected  geometrical 
enhancement  from  sheetlike  irregularities.  The  largest  scintillation 
enhancements,  however,  were  observed  on  GPS  8  after  1100  UT  where  the 
geometrical  enhancement  factor  is  minimized.  Coincident  Chatanika  radar 
data  have  revealed  a  persistent  enhanced  F-region  ionization  "wall"  that 
is  possibly  the  source  region  of  the  dominant  enhancement. 
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Table  1 

GPS  INTENSITY  SCINTILLATION  SUMMARY 


Sa tellite 

s4  >  0.1 

Maximum 

No. 

(7.) 

s4 

309  I 

8 

4 

310 

8 

4 

5 

311 

8 

4 

5 

313 

8 

4 

5 

314 

8 

4 

5 

315 

8 

4 

5 

316 

8 

4 

5 

317 

8 

5 

(College  K) 


5 

5,  4 

4  (5  after  midnight) 


6  (isolated  substorm) 
4 


3  (8  after  midnight) 


Computation  of  Phase  Turbulence  Strength 


We  showed  in  Section  I  that  the  beam  distortion  in  a  DFA  depends 
critically  on  the  phase  structure  constant  as  defined  by  Eq.  (22).  T 
scintillation  model  that  we  have  used  extensively  is  based  on  a  three 
dimensional  in  situ  spectrum  of  the  form 


where  q  is  a  quadratic  form  on  the  vector  elements  k^,  k  ,  and  k^.  If 
Eq .  (24)  is  integrated  over  all  spatial  wave  numbers,  the  rms  electron 
density  can  be  computed  as 


<AN 

e 


2 


> 


q 


-  (2v 

o 


2)  r(v  -  i ) 
8*3/2  r(\)  +  \) 


(25) 


Thus,  if  a  value  of  q  is  specified,  and  AN^  are  interchangeable. 
As  a  practical  matter,  the  spatial  homogeneity  required  for  a  meaningful 
spectral  definition  is  generally  violated  before  an  outer-scale  cutoff 
can  be  identified.  Nonetheless,  power-law  segments  such  as 


csq 


- ( 2v  +  1  ) 


q  .  I  q  <  q 

;mn  max 


(26) 


can  be  deduced  from  measurable  one-dimensional  spectra. 

For  example,  an  in  situ  probe  measures  a  one-dimensional  spectrum  of 
the  form 


where 


and 


*(1) 

$ 

AN 


(q) 


cgr(v  -  \) 

4s2r (v  +  \) 


=  2v  -  i 


(27) 


(28) 


(29) 


The  actual  value  of  q  depends  on  the  scan  rate  and  the  scan  direction. 
This  relation  was  used  by  Livingston  et  al.  (1981)  to  compare  in  situ 
AE-E  satellite  data  with  phase  scintillation  data  from  the  Wideband 
sate  1 1 ite  . 

For  phase  scintillation  data,  we  measure  a  temporal  phase  spectrum 
of  the  form 


<p(f)  =  Tf-p 


(30) 
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where 


T  =  C  G 


Vs"  r(v ) 


2v  -1 


p  (2s  )2v  +  1  r(v  +  '*)  eff 


(31  ) 


=  2v 


(32  ) 


and 


2  2 

C  =  r  X  £  C0 
P  e  P  S 


(33) 


As  discussed  in  Rino  (1982b),  Eq.  (22)  can  be  approximated  by  the  simpler 


4.5v  +  5.5) 


1.0  <  v  <  2.0 


(34) 


Hence,  once  and  v  are  determined,  can  be  computed. 

The  major  uncertainty  in  determining  from  T  is  the  effective  scan 
velocity  v  which  depends  on  both  the  irregularity  anisotropy  and  the 

apparent  irregularity  drift.  Because  the  auroral-zone  field  lines  are 
nearly  vertical,  for  rodlike  irregularities  v  is  simply  the  magnitude 
of  the  component  of  the  apparent  drift  velocity  perpendicular  to  the 
line  of  sight.  For  sheetlike  irregularities,  however,  the  east-west 
component  of  the  apparent  velocity  is  reduced  by  the  east-west  enhance¬ 
ment  factor.  The  apparent  velocity  itself  results  from  the  true  irregu¬ 
larity  drift  and  the  motion  of  the  ionospheric  penetration  induced  by 
the  signal  source  or  receiver.  For  the  GPS  observations,  the  source 
notion  is  negligible  because  of  the  12-hour  orbits. 

To  illustrate  the  interplay  between  anisotropy  and  the  temporal  to 

spatial  conversion  factor  v  rr,  a  plot  of  v  for  a  100-mps  westward 

eff  eff  r 

drift  is  shown  in  Figure  6.  The  two  sets  of  curves  are  for  5:1  rods  and 
5:5:1  sheets.  For  sheets,  the  temporal-to-spatial  conversion  factor  is 
fairly  uniform  and  the  true  drift  is  recovered  directly.  For  rods,  there 
is  much  more  variation.  With  a  single  point  measurement,  there  is  no  way 
to  resolve  such  ambiguities. 
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FIGURE  6  EFFECTIVE  VELOCITY  FOR  TEMPORAL  TO  SPATIAL  CONVERSION 


To  illustrate  the  general  characteristics  of  the  phase  scintillation 
data,  in  Figure  7  wo  show  the  phase  data  for  GPS  5  on  Day  310,  after  a 
300-s  detrend  filter  was  applied.  Spectral  analysis  was  then  performed 
on  iS-min  overlapping  segments.  Two  representative  spectra  are  shown  in 
Figure  7.  Tt  is  important  to  note  that  the  aggregate  disturbance  for  the 
earlier  spectrum  is  larger  than  for  the  second  disturbance,  even  though 
the  latter  shows  a  larger  peak-phase  disturbance.  Peak-to-peak  phase 
excursions  in  excess  of  20  radians  are  associated  with  all  the  moderately 
disturbed  data  segments.  That  fact  alone,  however,  gives  n '  iseful  in¬ 
formation  about  disruptive  system  effects. 

The  spectral  data  have  been  smoothed  and  summarized  bv  the  power-law 
parameters  T  and  p  obtained  from  a  log- linear  least-squares  fit  over  the 
temporal  frequency  range  0.01  to  0.03  Hz.  On  average,  the  spectral  slopes 
are  steeper  than  those  measured  using  Wideband  satellite  data.  This  is 
very  likely  from  temporal  variations  within  the  comparatively  long 
(30-min)  data  segments.  The  GPS  data  show  a  systematic  decrease  in  the 
spectral  slope  with  increasing  perturbation  strength  as  has  been  consis¬ 
tently  observed  in  the  Wideband  satellite  data.  This  is  illustrated  in 
the  p  vs.  log^T  plot  shown  in  Figure  8,  which  summarizes  the  GPS  5  data 
for  four  evenings.  The  straight  line  was  hand  drawn.  Similar  plots  for 
Wideband  satellite  data  show  the  same  trend,  but  displaced  to  a  median 
value  of  “V  -2.5  rather  than  ^  -3.0. 
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FIGURE  7  GPS  PHASE  SPECTRA  (DAY  310) 


For  equatorial  data  there  is  a  consistent  "V  5  : 1  ratio  between  the 
effective  scan  velocities  for  the  low-orbiting  Wideband  satellite  and 
GPS.  In  the  Wideband  satellite  phase  scintillation  data,  the  temporal 
frequencies  that  are  least  affected  by  diffraction  generally  correspond 
to  scale  sizes  greater  than  one  kilometer.  Thus ,  the  Wideband  satellite 
equatorial  phase  data  are  heavily  weighted  toward  the  "intermediate- 
scale  regime"  where  the  slopes  are  shallow  and  vary  systematically 
[Livingston  et  al.,  1981].  The  corresponding  range  of  scale  sizes  for 
the  GPS  equatorial  data  possibly  lies  below  one  kilometer  where  the 
spectrum  shows  a  distinctly  different  behavior  [Rino  et  al.,  1981]. 

In  the  auroral  zone,  the  in  situ  drifts  are  substantially  larger 
and  more  variable  than  the  equatorial  drifts.  Thus,  we  do  not  tiave  a 
sharp  separation  of  the  scale-size  regimes  sampled  by  Wideband  and  by 
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FIGURE  8  SPECTRAL  INDEX  VARIATION 

GPS;  moreover,  because  the  auroral-zone  is  much  richer  in  plasma  source 
mechanisms,  irregularity  drivers  and  coupling  regions,  the  spectral 
characteristics  themselves  are  likely  to  be  more  varied. 

C.  Single  Baseline  Spaced- Receiver  Measurements 

To  complete  the  analysis  of  DFA  beam  smearing,  we  need  to  know  the 
anisotropy  of  the  diffraction  pattern.  To  resolve  the  anisotropy,  unam¬ 
biguously,  multireceiver  measurements  are  required  [Ring  and  Livingston, 
1982],  On  the  other  hand,  a  single  spaced-recei ver  measurement  is  the 
simplest  DFA  configuration  that  retains  all  the  essential  elements  of 
the  system.  In  the  plane  of  the  DFA,  the  correlation  function  of  the 
field  admits  a  representation  of  the  form 

R[  f  ( Ap*'  -  v'  6  t )  ]  (35) 


where 


f2(C)  =  \c.X2  -  B'C  (  +  A'C21/[A'C'  -B'2/4]  .  (36) 

Lx  x  y  yj 
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As  discussed  in  Rino  and  Livingston  [1982],  the  anisotropy  of  the  signal 
intensity  and  the  signal  phase  is  essentially  the  same  as  the  complex 
signal  itself;  however,  because  of  the  small  perturbations  levels  at 
L  band,  we  shall  use  only  signal  phase. 

By  operating  two  independent  dual-frequency  GPS  receiver  systems, 

we  were  able  to  make  d isplaced-point  measurements  of  the  signal  phase 

from  a  common  source.  Because  the  transmission  time  is  modulated  onto 

the  data  stream  and  because  the  pseudo-range  delay  differences  are  at 

most  many  tens  of  nanoseconds  even  without  bias  removal ,  the  phase 

measurements  can  be  coherently  processed  for  phase  comparisons.  Indeed, 

this  is  the  phase  error  that  a  100-m  linear  DFA  would  encounter  if  it 

were  "focused"  on  the  GPS  satellite.  Because  the  error  is  random,  the 

degree  of  phase  decorrelation  between  the  two  receivers  is  the  quantity 

of  interest;  call  it  R,(Ap). 

0 

We  have  shown,  however,  that  the  phase  structure  function,  Eq .  (20), 
is  the  quantity  that  characterizes  the  average  DFA  beam  distortion.  The 
formal  interrelation  is 

Dfi0(Ap)  =  2[R^(0)  -  R^fAp)]  ,  (37) 

provided  that  diffraction  effects  in  the  signal  phase  are  negligible. 

The  effect  of  Eq.  (37)  is  to  suppress  the  contribution  of  large-scale 

— ♦ 

structures.  Treatments  that  proceed  directly  from  R^(Ap )  invariably 
extract  these  low  frequency  components  [Elachi  and  Evans,  1977], 

If  a  coordinate  system  with  one  axis  along  the  receiver  baseline  is 
used,  Eq.  (35)  can  be  rewritten  as 

R[  (£,  -  v  6 1 )  /£  ]  (38) 

b  a  c 

where  £,  is  the  baseline  length,  v  is  the  apparent  velocity,  and  £  is 
d  a  c 

the  apparent  correlation  scale.  The  problem  with  single  baseline  measure- 

ments  is  that  £  and  v  are  not  simply  related  to  the  minimum  spatial 
c  a 

coherence  scale  and  the  effective  or  true  pattern  drifts. 
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Unfortunately,  the  data  acquisition  system  (which  was  rapidly  con¬ 
structed  just  before  the  experiment)  for  the  second  receiver,  did  not  perform 
well,  and  we  obtained  only  a  few  coincident  data  sets  during  interesting  scin 
dilation  events.  These  are  summarized  in  Table  2.  The  first  two  columns 
list  the  maximum  value  of  the  cross-correlation  function  and  its  delay  offset 
respectively.  The  remaining  columns  summarize  the  phase  spectral  character¬ 
istics  at  the  individual  receivers.  The  fact  that  the  spectral  character¬ 
istics  are  identical  to  within  the  statistical  accuracy  of  the  measurements 
is  a  good  indication  that  all  the  instrumental  effects  have  been  removed. 

The  actual  phase  records  and  the  measured  cross  correlation  functions 

for  days  318  and  320  are  shown  in  Figures  9  and  10,  respectively.  The 

data  in  Figure  9  were  taken  during  active  positive  hay  conditions  as 

recorded  on  the  College,  Alaska  magnetometer.  For  T  'V  1  and  £,  'V  100  m, 

a  b 

v  is  approaching  1  km/s.  Unlike  the  effective  velocity  for  converting 

temporal  scales  to  spatial  scales  (v^^.),  the  apparent  pattern  drift  will 

always  be  larger  than  the  true  drift.  Thus,  large  values  of  v^  are  more 

indicative  of  enhanced  correlation  scales  along  the  baseline  that  increases 

in  the  pattern  drift.  The  t  values  approaching  0.9,  however,  correspond 

a 

to  several  hundred  meter-per-second  east-west  drifts  which  are  consistent 
with  Chatanika  radar  drift  measurements  made  under  similar  conditions. 

The  data  shown  in  Figure  10  are  more  interesting  in  that  they  show 
an  abrupt  change  in  activity  around  928  UT.  An  abrupt  negative  bay  on¬ 
set  was  observed  on  the  College  magnetometer  just  before  900  UT.  Because 
the  GPS  5  satellite  was  almost  10°  to  the  west  of  Poker  Flat,  this  event 
probably  marks  the  passage  of  the  westward  expansion  of  the  activity  to 

the  ionospheric  penetration  point.  The  eastward  drift  component  (v  <  0) 

a 

of  A.  600  mps  is  reasonable  for  an  active  eastward  electrojet. 

Thus  the  GPS  dua 1 -recei ver  measurements  seem  to  produce  credible 
measurements  that  can  be  exploited  further.  Insofar  as  DFA  performance 
is  concerned,  for  the  perturbation  levels  experienced  here,  there  was 
essentially  no  phase  decorrelation  on  a  1 00-m  baseline.  It  follows  that 
small  DFA  systems  operating  at  L  band  or  higher  would  experience  no  sig¬ 
nificant  beam  degradation.  This  will  be  shown  quantatively  in  Section  5 
in  which  all  the  results  are  discussed  and  conclusions  drawn. 


Table  2 

GPS  SPACED- RECEIVER  SUMMARY 


Day 

m 

Maximum 

Correlation 

Delay 

(s) 

ti 

(dB) 

T2 

(dB) 

P1 

P2 

M 

0741  to  0807 

0.12 

ra 

-  2.6 

-  2.7 

0754  to  0819 

KEffl 

0.94 

Wmm 

-  2.7 

-  2.6 

0807  to  0832 

0.92 

0.80 

-  48 

-  2.5 

-  2.8 

319 

1451  to  1517 

■I1M 

-  0.13 

-  33 

-  32 

-  2.2 

-  2.1 

(GPS-6) 

1503  to  1528 

-  0.05 

-  26 

-  28 

-  1.8 

-  1.9 

1516  to  1541 

ks 

-  0.02 

-  32 

-- 

-  2.3 

— 

320 

0905  to  0930 

mm 

-  55 

-  60 

-  3.6 

-  3.9 

(GPS-5) 

0912  to  0937 

■m 

Bail 

-  42 

-  45 

-  3.1 

-  3.2 

1.5 


FIGURE  9  GPS  SPACED- RECEIVER  MEASUREMENTS  (DAY  318) 
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FIGURE  10  GPS  SPACED-RECEIVER  MEASUREMENTS  (DAY  320) 
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V  MULTIPLE  BASELINE  SPACED- RECEIVER  ANALYSIS 


A.  Introduction 

In  the  conclusion  of  Section  II,  we  pointed  out  that  an  essentia] 
element  for  calculating  DFA  average  beam  distortion  is  irregularity 
anisotropy  because  the  severity  of  a  propagation  disturbance  is  highly 
aspect  sensitive.  The  only  way  to  unambiguously  establish  this  anisotropy 
at  the  spatial  scales  important  for  DFA  operation  is  through  the  use  of 
multiple-baseline  spaced  receivers.  These  same  measurements  also  provide 
the  irregularity  drift  information  required  in  the  calculation  of  temporal 
perturbation  strength.  Thus,  if  phase  scintillation  is  measured  concur¬ 
rently,  as  it  was  for  the  data  discussed  here,  the  final  DFA  analysis 
element,  turbulence  strength,  can  also  be  determined. 

The  application  of  the  spaced-receiver  method  to  measure  anisotropy 
is  straightforward.  In  situ  irregularities  can  generally  be  characterized 
in  terms  of  a  three-dimensional  spatial  autocorrelation  function;  the 
radio-wave  pattern  produced  at  the  receiver  plane  by  illuminating  those 
irregularities  will,  therefore,  have  a  prescribed  two-dimensional  form. 

For  example,  if  the  in  situ  irregularities  were  prolate  spheroids,  they 
will  produce  a  diffraction  pattern  with  elliptical  contours  of  constant 
correlation.  By  making  multiple,  simultaneous  spatial  samples  of  that 
diffraction  pattern  at  the  receiver  plane,  correlation  techniques  can  be 
used  to  characterize  its  shape  and  orientation.  By  reversing  the  process, 
these  parameters  can  be  used,  in  conjunction  with  a  propagation  model,  to 
determine  the  in  situ  structure. 

Similarly,  drift  measurements  using  multiple  receivers  are  simple  in 
concept.  A  frozen  diffraction  pattern  moving  across  the  receiver  plane 
will  produce  identical  signals  displaced  in  time.  In  the  simplest  case 
of  an  isotropic  correlation  pattern,  dividing  the  baseline  separation  by 
the  time  delay  yields  the  drift  component  along  that  baseline. 


t 
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In  reality,  the  contributions  of  anisotropy  and  drift  to  the  spatial 
correlation  pattern  are  not  directly  separable,  as  discussed  in  Section  III. 
The  interplay  between  anisotropy  and  drifts  for  one  measurement  geometry 
is  shown  in  Figure  11.  Two  differently  shaped  irregularities  are  illumi¬ 
nated  at  an  angle  relatively  close  to  the  direction  of  magnetic  field 
(i.e.,  small  Briggs-  Parkin  angle ).  On  the  right,  the  irregularity  is 
field-aligned  and  rodlike.  Projection  of  the  irregularity  pattern  onto 
the  receiving  plane  produces  a  correlation  surface  with  anisotropy  lower 
than  that  in  situ,  and,  in  this  case,  an  alignment  near  the  magnetic 
meridian.  On  the  left,  the  irregularity  has  a  second  preferred  axis  along 
the  local  L  shell,  resulting  in  a  sheetlike  structure.  Because  of  this 
magnetic  east-west  extension,  the  ground  pattern  is  more  highly  anisotropic 
than,  and  extended  in  a  direction  perpendicular  to,  that  for  rods. 

This  example  shows  how,  from  a  correlation  surface  characterized  in 
terms  of  axial  ratio  and  orientation,  the  two  in  situ  irregularity  types 
can  be  readily  distinguished.  Equally  important,  the  example  also  under¬ 
scores  the  need  to  account  for  pattern  anisotropy  before  drift  information 
can  be  extracted  from  spaced-receiver  data.  The  difficulty  is  that  a  high 
degree  of  correlation  along  a  particular  baseline  will  be  interpreted  as 
a  large  pattern  drift  if  one  simply  divides  the  correlation  delay  into 
the  baseline  length.  For  rods,  the  signals  at  Antennas  1  and  2  are  highly 
correlated  with  peaks  only  slightly  displaced  in  time.  Interpretation 
of  that  small  temporal  delay  as  drift  produced,  without  proper  correction 
for  pattern  anisotropy,  would  suggest  a  high  north-south  drift.  Similarly, 
the  sheet-produced  pattern  would  indicate  high  east-west  velocities,  which 
can  become  infinite  if  the  pattern  and  baseline  are  exactly  aligned,  and 
would  be  observed  even  if  the  irregularities  were  not  moving.  True  drift 
velocity  information  can  be  extracted  only  if  the  spatial  pattern  aniso¬ 
tropy  is  accounted  for. 

Tn  the  results  presented  in  the  next  section,  we  have  used  a  corre¬ 
lation  method  from  which  pattern  anisotropy  and  true  pattern  drift  can 
be  simultaneously  extracted.  The  technique  based  upon  the  work  of 
Armstrong  and  Coles  fl972]  and  is  described  in  Rino  and  Livingston  fl982]. 
The  method  uses  the  intersections  and  peaks  of  all  possible  cross-  and 

auto-correlation  fvinctions  derived  from  three  or  more  antennas.  An 
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over-determined  act  of  linnar  equations  results  that  In  solved  for 
anisotropy  n nci  true*  pattern  drift  by  least-squares  methods.  Tho  observed 
values  arc  thon  compared  with  modol  predictions  to  establish  the  anlso- 
tropy  of  the  scintillation-producing  structures.  Tho  technique  can  bo 
appllod  to  signal  intensity  or  to  signal  phase.  The  effoctn  of  diffrac¬ 
tion  alter  the  overall  size  of  tho  spatial  correlation  function*  of  tho 
signal  at  tho  receiver,  but  they  do  not  alter  Its  anisotropy. 

Generally  spooking,  the  parameters  that  describe  the  signal  spatial 
correlation  at  the  receiver  arc  unique  for  a  particular  generic  type  of 
irregularity,  in  our  analysis  we  consider  three  such  typos.  Fiold- 
allgned  rods  con  be  described  in  terms  of  an  axial  ratio,  a,  which  is  the 
ratio  of  tho  averago  along-fleld  dimension  to  the  cross-field  dimension, 
and  are  designated  a:1:l.  When  tho  irregularities  have  a  socond  preferred 
axis  perpendicular  to  the  field  and  along  the  local  L-shcll,  they  arc 
referred  to  as  sheets  and  are  designated  a:a:l.  Similarly  orlontcd,  but 
nonsyrnmctrlcol  Irregularities,  (nib:  I  where  b  <  a)  are  referred  to  as 
wings . 

For  application  to  the  DFA  analysis,  we  have  used  data  from  two  very 
different,  but  complementary,  signal  sources.  The  first  of  these  is  the 
Wideband  satellite  f Fromouw  et  a  1 . ,  1978],  Most  of  the  data  from  that 
experiment  has  boon  already  published,  and  we  will  present  only  summary 
examples  here.  Its  sun-synchronous,  106-min  orbit  provided  latitudinal 
scan  snapshots  of  the  auroral  zone,  within  a  few  hours  of  local  midnight, 
over  soveral  seasons  of  low-to-modorato  sunspot  conditions.  Tho  other 
signal  sources  that  wo  have  utilized  are  a  sorics  of  UHF  beacon  satellites 
in  highly  occontrlc  orbits  and  can  be  observed  continuously  for  a  few  hours 
at  northern  high  latitudes.  These  satellites  have  boon  observed  from 
stations  in  the  auroral,  zone  and  polar  cap  routinely  since  1979,  although 
we  will  concentrate  ori  data  collected  during  recent  coordinated  measure¬ 
ment  campaigns.  The  nature  of  the  data  from  thoso  latter  measurements  is 
quite  different  from  the  0P8  or  Wideband  satellite  data.  For  example, 
tho  penetration  point  moves  very  little,  so  that  the  geometrical 
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propagation  factor  G  (Eq.  31)  remains  nearly  constant  and  the  contribution 
of  satellite  motion  to  the  apparent  irregularity  drift  is  small,  being 
confined  primarily  to  the  radial  direction.  The  temporal  behavior  of 
anisotropy  and  drift  at  a  nearly  constsnt  latitude  can,  therefore,  he 
observed  daily  over  a  several  hour  window. 

We  have  presented  our  examples  in  two  sections.  Section  V-R  discusses 
the  auroral  zone  and  concentrates  on  anisotropy  in  the  evening  and  midnight 
sectors  and  a  comparison  between  our  drift  measurements  and  those  made  by 
the  Chatanika  incoherent-scatter  radar.  Turbulence  strength  calculations 
for  strong  auroral  zone  scintillation  are  also  presented.  Section  V-C 
presents  examples  of  anisotropy,  drifts  and  turbulence  strength,  measured 
in  the  central  polar  cap,  which  provide  guidelines  for  DFA  applications 
at  those  latitudes. 

B.  Measurements 

1 .  Auroral-Zone  Irregularity  Anisotropy  and  Drifts 

Kilometer-scale,  scintillation-producing  irregularities  in  the 
auroral  zone  have  been  observed  to  have  a  highly  variable  anisotropy.  At 
these  scales,  scintillation  is  enhanced  over  propagation  paths  lying  along 
any  extended  axis  of  an  irregularity  [Rino  and  Fremouw,  1*377],  a  phenomenon 
that  has  been  used  as  a  diagnostic  in  comparing  measured  and  predicted 
scintillation  strength  as  a  function  of  propagation  geometry.  In  the 
auroral  zone,  localized  enhancements  are  regularly  seen  not  only  whore 
the  propagation  is  along  the  magnetic  field,  as  produced  by  rodlike  struc¬ 
tures,  but  also  when  the  path  lies  in  the  local  L  shell.  Martin  and  Aarons 
[1977]  explained  this  occurrence  by  postulating  a  second  preferred  irregu¬ 
larity  axis,  i.e.,  L-shel 1-al igned  sheets.  Rino  anti  Matthews  [1980]  and 
Basu  et  al.  [1981]  arrived  at  a  similar  conclusion  using  scintillation 
enhancements  from  the  extensive  Wideband  satellite  data  base.  Observa¬ 
tions  from  both  the  Pacific  and  Atlantic  sector  auroral  zone  were  included 
in  these  latter  studies.  Unfortunately,  the  method  depends  upon  scintil¬ 
lation  enhancements  being  totally  attributable  to  geometry.  Accordingly, 
the  results  can  be  ambiguous  if  irregularity  strength  is  not  constant  with 
latitude  and  time. 
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Moore ro ft  and  Arima  [1972]  were  the  first  to  determine  anisotropy 
from  spaced -rec e i ver  data  in  the  auroral  zone.  They  found  evidence  of 
sheetlike  structures,  hut  with  meridional  alignment.  More  recently, 

Rino  and  Livingston  [1982]  present  evidence  of  sheetlike  and  winglike 
structures,  derived  from  Wideband  spaced-recei ver  data.  This  confirmed 
the  occurrence  of  sheets  demonstrated  earlier  by  Rino  et  al.  [1978], 
based  on  spaced-recei ver  analysis  of  data  recorded  at  Poker  Flat,  Alaska. 

A  more  complete  study  of  anisotropy  was  recently  made  by 
Livingston  et  al.  [1982],  also  based  upon  Wideband  data.  They  established 
that  auroral-zone  anisotropy  depends  upon  both  latitude  and  local  time 
and  appears  to  be  convection  dominated.  Figure  12  is  a  reproduction  of 
one  of  their  summary  plots,  which  corresponds  to  moderate  magnetic  and 
sunspot  number  conditions.  In  the  pre-  and  post-midnight  sectors,  the 
nor th-to-south  latitudinal  progression  is  from  nearly  isotropic  irregular¬ 
ities  to  rods  with  increasing  elongation.  These  give  way  to  sheets  in  a 
narrow  latitudinal  sector  corresponding  to  the  equatorward  edge  of  the 
aurora.  dust  before  local  midnight,  the  anisotropy  of  the  rods  lessens, 
and  the  sheets  become  wings.  Livingston  et  al.  [1982]  interpret  this 
pattern  by  the  Harang  discontinuity  f  Heppner ,  1972].  To  indicate  this 
process,  the  large-scale  flow  pattern  expected  for  this  level  of  magnetic 
activity  has  been  superimposed.  We  emphasize  that  this  is  a  model  pattern. 

Figure  12  suggests  that  the  sheets  and  wings  are  most  prevalent 
within  an  hour  or  two  of  local  midnight.  Unfortunately,  the  Wideband 
orbit  was  such  that  the  presence  of  sheets  in  the  evening  or  morning 
sectors  can  be  neither  confirmed  nor  denied.  (It  has  been  established, 
however,  that  during  the  hours  around  local  noon  no  sheets  or  wings  are 
observed .  ) 

The  question  of  whether  sheets  occur  in  the  evening  or  morning 
auroral  zone  is  important  both  to  understanding  of  irregularity  generation 
mechanisms  and  the  UFA  pertinent  morphology.  The  data  from  the  high- 
eccentricity-orbit  signal  sources  previously  mentioned,  helps  fill  this 
local  time  gap.  For  the  November  1981  Chatanika  radar  campaign,  between 
1800  LT  and  2100  LT,  the  geometry  to  those  sources  was  favorable  for 
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detailed  anisotropy  and  drift  measurements .  During  the  more  recent, 
February-to-March  1982  campaign,  also  at  Chatanlka,  the  optimum  geometry 
occurred  In  the  local  morning  sector. 

Our  discussion  here  will  concentrate  on  the  first  of  these  cam¬ 
paigns,  during  which  the  GPS  data  already  discussed  was  collected.  This 
Includes  data  that  directly  overlap  the  Wideband  data  of  Figure  12,  so 
that  the  temporal  variation  of  anisotropy  and  drift  on  one  evening  can 
be  directly  compared  with  the  long-term  morphological  pattern.  We  have 
also  Included  measurements  that  were  coordinated  with  Incoherent-scatter 
radar  observations.  This  provides  an  opportunity  to  (1)  establish  the 
Irregularity  anisotropy  at  Intermediate  spatial  scales  and  associate 
changes  In  the  structure  with  large-scale  plasma  features  observed  by 
the  radar,  and  (2)  compare  radar-measured  plasma  drifts  with  those  of 
the  Intermediate  scale,  scintillation-producing  structures. 

We  will  first  consider  the  midnight  sector  beacon  measurements. 
Figure  13  shows  the  correlation  pattern  axial  ratio  and  orientation 
measured  on  Day  315,  1981,  for  a  low-elevation,  northward-looking  geometry. 
Superimposed  on  the  measurements  are  the  model  predictions  for  5:1:1  rods 
and  5:5:1  sheets.  For  this  geometry  it  is  impossible  to  distinguish  the 
two  irregularity  types  from  pattern  orientation  alone.  The  axial  ratio 
data,  however,  do  provide  information  about  the  in  situ  structure. 

Sheets,  illuminated  broadside  produce  a  near-isotropic  correlation  pattern, 
while  rods  have  an  extended  shadow.  The  measured  axial  ratio  data  is 
systematically  higher  than  that  for  sheets,  but  lower  than  that  for  5:1:1 
rods,  indicating  an  overall  in  situ  anisotropy  of  about  3:1:1. 

Because  they  overlap  in  latitude  and  local  time,  we  can  compare 
the  anisotropy  pattern  in  detail  with  the  Wideband  measurements.  The  in 
situ  anisotropies  implied  from  the  model-measurement  comparisons  have  been 
mapped  in  Figure  14  in  the  same  format  as  the  Wideband  data.  To  facili¬ 
tate  the  comparison  with  Figure  12  and  to  indicate  the  convection  process 
that  we  believe  controls  the  anisotropy,  the  expected  large-scale  flow 
pattern  is  superimposed  on  Figure  14.  Both  Figures  12  and  14  correspond 
to  moderate  magnetic  conditions  (k  *  3  to  4)  and  should  therefore  be 
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FIGURE  13  MEASURED  AND  MODEL  PATTERN  ANISOTROPY 


directly  comparable  in  terms  of  the  latitudinal  and  local  magnetic  time 
dependence  of  anisotropy.  We  do,  in  fact,  see  a  consistent  behavior  be¬ 
tween  the  two.  A  local-time  cut  through  Figure  12  at  an  invariant  lati¬ 
tude  of  70°  shows  a  pattern  of  rods  ranging  from  3:1:1  through  5:1:1  with 
short  peaks  of  8:1:1  elongation  near  2220  and  2350  MLT.  This  is  very 
similar  to  the  anisotropy  pattern  of  Figure  14  where  peaks  of  5:1:1  rods 
near  2230  and  0000  MLT  are  observed.  The  minimum  in  irregularity  elongation 
is  not  quite  as  distinct  in  the  Day  315  data,  nor  is  the  generally  higher 
anisotropy  post-  versus  pre-midnight.  Those  differences,  and  the  generally 
lower  values  of  anisotropy,  are  most  likely  statistical.  On  the  other 
hand,  it  m.iv  indicate  that  spatial  scale  sire  depends  slightly  on  anisotropy. 
For  Wideband,  the  scales  sampled  are  in  the  1-  to  2-km  range,  as  dictated 
by  the  phase  detrend  cutoff.  The  beacon  data  corresponds  more  nearly  to 
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0.5  km,  which  is  near  the  signal  intensity  Fresnel  radius  cutoff  for  the 
observation  frequency  and  geometry. 

In  Livingston  et  al .  [ 1982]  it  was  postulated  from  the  anisotropy 
pattern  that  the  irregularity  structure  was  controlled  by  convection.  With 
the  UHF  polar-beacon  data,  where  velocities  are  simultaneously  extracted, 
it  is  possible  to  test  this  hypothesis  directly.  The  drift  data  are  plotted 
in  Figure  15.  We  emphasize  that  these  are  true  (not  apparent)  drift  veloci¬ 
ties  from  which  the  pattern  anisotropy  has  been  removed.  Superimposed 
on  Figure  15  is  the  same  Harang  discontinuity  pattern  used  in  Figures  ]3 
and  14.  As  can  be  seen,  the  data  in  Figure  15  match  the  predicted  large- 
scale  flow  pattern  remarkably  well.  Before  2300  CMLT,  the  drift  is  about 
250  m/s  eastward  and  switches  to  about  300  m/s  southward  by  local  midnight. 

To  demonstrate  auroral  zone  anisotropy  and  drifts  in  a  local 
time  sector  which  Wideband  did  not  sample,  we  now  turn  to  data  collected 
from  a  different  polar  beacon  early  in  the  evening  of  11  November  1981. 

The  F-region  penetration  latitude  for  these  data  is  near  65°  invariant  or 
near  the  poleward  edge  of  the  sheet-wing  region  in  Figure  13,  although  as 
we  have  said,  the  two  experiments  do  not  overlap  in  local  time. 

Figure  16  shows  the  measured  correlation  pattern  anisotropy  and 
orientation  for  this  evening,  and  the  superimposed  model  predictions  for 
5:1:1  rods  and  5:5:1  sheets.  Unlike  the  lower  elevation  geometry  of 
Figure  13,  the  rods  and  sheets  produce  patterns  that  are  significantly 
different  in  orientation.  Note  that  extending  the  rod  and  sheet  dimen¬ 
sions  beyond  that  used  here  will  not  alter  the  correlation  surface  orien¬ 
tation,  although  it  will  alter  the  observed  pattern  axial  ratio.  Figure 
16  indicates  that  during  the  observation  period  the  in  situ  irregularities 
are  rodlike.  At  no  time  does  the  orientation  deviate  systematically  from 
the  rod  pattern  towards  that  for  sheets.  The  axial  ratio  does  vary,  how¬ 
ever,  corresponding  to  in  situ  rod  extensions  from  a  minimum  of  5:1:1  to 
a  maximum  of  15:1:1.  Incidentally,  these  changes  seem  to  be  associated 
with  the  passage  of  macroscale  plasma  enhancements  seen  in  the  radar  data, 
although  a  detailed  correlation  has  yet  to  be  made. 

Because  these  and  the  Wideband  data  cannot  he  directly  compared, 
it  is  difficult  to  say  whether  the  absence  of  sheets  in  the  polar  beacon 
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FIGURE  15  MEASURED  IRREGULARITY  DRIFTS 
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FIGURE  16  MEASURED  AND  MODEL  PATTERN  ANISOTROPY 

data  can  be  attributed  to  a  local  time  or  latitudinal  dependence.  Wo  note, 
however,  that  the  in  situ  anisotropy  is  higher  for  this  near-overhead 
geometry  than  to  the  north  later  in  the  evening.  In  this  sense,  the  lati¬ 
tudinal  behavior  of  the  degree  of  rod  elongation  is  consistent  between  the 
evening  sector  and  the  midnight  sector  data  of  Figure  12. 

An  advantage  of  the  near-overhead  geometry  for  this  data  set  is 
that  the  simultaneously  extracted  irregularity  drifts  can  be  compared  with 
radar-measured  plasma  drifts.  The  spaced-recei ver  drift  data  is  shown  in 
Figure  17,  rotated  into  magnetic  coordinates.  These  are,  again,  corrected 
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FIGURE  17  MEASURED  IRREGULARITY  DRIFT  COMPARISON  WITH  CHATANIKA  RADAR 
AND  DE-B  SATELLITE  DRIFT  DATA 

for  pattern  anisotropy  and  peak  at  nearly  800  m/s  westward  near  1900  CMLT. 
The  corresponding  apparent  velocities,  that  is,  those  obtained  by  dividing 
the  correlation  peak  delays  into  the  baseline  length,  routinely  exceed 
2000  m/s.  Clearly,  interpretation  of  the  apparent  velocities  as  true 
drifts  leads  to  serious  error. 

The  Chatanika  radar  F-region  plasma  drifts  are  superimposed  on 
Figure  17.  Those  values  arc  those  measured  at  the  same  penetration  lati¬ 
tude  as  the  satellite  data  and  are  derived  by  combining  elevation  scan 
pairs:  one  in  the  meridian  and  one  to  the  west  at  the  same  magnetic 
local  time.  As  the  map  in  Figure  18  illustrates,  these  scans  bracket 
the  satellite  penetration  point  in  longitude;  the  duration  of  the  scans 
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FIGURE  18  MAP  OF  THE  MEASUREMENT  GEOMETRY  FOR  COMPARISON 
OF  SPACED-RECEIVER,  RADAR,  AND  IN-SITU  DRIFT 
MEASUREMENTS 

gives  a  single  velocity  about  every  25  min,  while  those  from  the  scintil¬ 
lation  measurements  are  about  every  1.5  min. 

The  north-south  drift  is  low  as  expected  for  the  local  evening 
sector.  The  two  measurements  match  well,  although  the  spaced-rece i ver 
system  provides  much  more  temporal  resolution.  The  radar  indicates  a 
strong  westward  drift,  which  eventually  exceeds  800  m/s  starting  about 
1700  CMLT,  and  which  corresponds  to  the  passage  of  a  500-y  positive  bay. 
The  spaced-receiver  drifts  show  a  very  similar  temporal  pattern,  although 
initially  the  westward  irregularity  drift  is  much  smaller  in  magnitude 
than  that  measured  by  the  radar.  The  biggest  difference  between  the  two 
measurements  occurs  near  1715  CMLT,  when  the  westward  irregularity  drift 
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1 8  nhout  60  percent  of  the  plasma  drift  magnitude;  by  1900  CMLT,  tbo  two 
measurements  match  very  well.  Two  oxplanatlonH  can  be  offered  for  the  early 
time  discrepancy.  The  validity  of  tho  radar  measurement  depends  upon 
homogeneity  across  the  several  degrees  of  longitude  separating  the  scans. 
This  condition  may  not  be  satisfied  during  the  Initial  passage  of  the  mag¬ 
netic  disturbance,  leaving  the  radar  measurements  questionable.  After 
about  1830  CMLT,  when  the  disturbance  has  passed  and  longitudinal  homo¬ 
geneity  Is  restored,  the  measurements  agree. 

Alternatively,  the  Irregularity  drifts  can  Include  the  effects 
of  local  polarization  fields  If  the  density  enhancement  within  the  Irregu¬ 
larity  Is  strong.  For  example,  If  the  north-south  polarization  field  pro- 

V 

duced  by  the  Irregularities  were  In  the  order  of  10  m  /m,  the  westward 
drift  of  the  Irregularities  at  our  600  m  spatial  scale  would  be  reduced 
by  about  150  m/s  from  that  cf  the  plasma.  Under  the  right  conditions, 
a  large-scale  density  enhancement  of  3  *  10^  could  alter  the  field  this 
much  (Tsunoda,  private  communication)  although  quantitative  estimates  of 
the  effect  would  be  difficult,  particularly  the  effects  of  lower  ionosphere 
conductivity.  For  either  mechanism,  the  scenario  agrees  well  with  the 
plasma  conditions  seen  by  the  radar.  Up  through  1715  CMLT,  the  iono¬ 
sphere  is  relatively  uniform.  Between  that  time  and  1900  CMLT,  a  number 
of  100-km  scale  density  enhancements  structured  in  latitude  and  longitude 
drift  through  the  sampled  region  and  then  disappear  to  the  south. 

A  second  velocity  comparison  can  be  made  with  these  data  using 
the  Dynamics  Explorer  satellite  (DE-B),  which  passed  the  penetration 
latitude  to  the  west  of  Chatanika  as  shown  in  Figure  18.  The  in  situ 
measurement  technique  is  described  by  Heells  et  al.  [1981].  During  the 
satellite  pass,  which  occurred  near  1805  CMLT,  we  can  assume  that  the 
horizontal  drift  component  is  purely  magnetic  westward  from  the  satellite 
orbital  inclination  and  the  dominance  of  the  westward  drift.  The  magnitude 
of  the  DE-B  measured  drift  is  about  600  m/s,  although  velocity  increases 
a  few  degrees  south  in  latitude.  This  single  comparison  point  is  indi¬ 
cated  in  Figure  17  and  matches  the  spaced  receiver  derived  value  almost 
exactly.  Because  the  in  situ  probe  responses  to  local  fields,  better 
agreement  with  the  spaced-receiver  measurement  than  with  the  radar  data 
could  be  expected. 


44 


In  conclusion,  the  auroral  zone  spaced-recei ver  measurements 
have  provided  anisotropy  and  drift  data  directly  applicable  to  the  DFA 
problem.  The  validity  of  the  drift  extraction  technique  has  been  proven 
by  comparison  with  radar  and  in-situ  data.  In  the  midnight  sector,  the 
anisotropy  variations  with  latitude  and  local  time  have  been  established 
using  Wideband  data,  and  verified  by  the  polar  beacon  observations.  In 
the  early  evening  sector,  the  temporal  behavior  of  anisotropy  in  the  central 
auroral  zone  has  been  measured.  The  agreement  between  our  drift  measure¬ 
ments  and  classical  flow  patterns,  and  the  apparent  convective  control  of 
anisotropy  suggests  that  the  application  of  a  global  model  for  DFA  analyses 
should  be  straightforward. 

2 .  Auroral-Zone  Turbulence  Strength  and  DFA  Limitations 

The  night  of  11  November  1981,  on  which  the  anisotropy  and 
velocity  observations  in  Figure  4  and  5  were  made,  was  also  a  period  of 
strong  intensity  and  phase  scintillation.  The  intensity  scintillation 
index  at  244  MHz  was  in  or  near  saturation  for  two  hours  near  midnight  and 
roughly  corresponds  to  the  maximum  scintillation  levels  observed  by  the 
Wideband  experiment  at  the  same  propagation  geometry.  If  coupled  with 
field-aligned  propagation  geometries  and  high  plasma  velocities,  this 
strong  irregularity  turbulence  may  represent  near-worst  case  condition 
for  DFA  application  in  the  auroral  zone. 

The  starting  point  for  the  calculation  of  turbulence  strength, 

Cg ,  (Eqs.  30  through  33)  is  the  measured  phase  spectral  density.  The 
strongest  scintillation  periods  on  II  November  correspond  to  average 
phase  spectral  strength,  T  =  -|8  dB,  and  spectral  index,  p  =  2.9.  In 
addition  to  these  measured  parameters,  we  have  assumed  (1)  rodlike  irreg¬ 
ularity  anisotropies  of  between  5:1:1  and  10:1:1,  typical  of  those 
measured,  (2)  a  plasma  drift  velocity  of  800  m/s,  which  is  near  the 
maximum  observed  during  the  period,  and  (3)  a  scattering  layer  thickness 

of  100  km.  For  the  propagation  geometry  of  the  observations,  this  results 

20 

in  a  value  of  C  =  10  It  is  interesting  to  note  that  this  level  is 

s 

near  the  median  value  observed  at  the  magnetic  equator  during  moderate 
sunspot  number,  disturbed  evening  conditions  [ Li vi ngs ton  e t  a  1 . ,  1981]. 

This  is  consistent  with  the  Wideband  satellite  observations,  which 


showed  consistently  weaker  scintillation  ir.  the  auroral  zone  than  at  the 
equator,  particularly  as  the  sunspot  activity  increased. 


The  degree  of  DFA  beam  distortion  can  be  estimated  by  calculat¬ 
ing  the  spatial  decorrelation  contour  described  by  Eq  (21).  We  consider 
two  propagation  situations.  The  first  of  these  is  for  the  observed 
average  axial  ratio  of  5:1:1  rods  and  the  actual  observation  geometry. 

The  second  case  is  for  field-aligned  propagation  and  the  peak  observed 
axial  ratio  of  10:1:1.  This  second  situation  is  worse-case  in  two 
regards.  The  propagation  factor  G  peaks  at  field  alignment  so  that  the 
perturbation  strength  maximizes.  In  addition,  the  spatial  correlation 
pattern  is  nearly  isotropic  so  that  decorrelation  is  equally  rapid  in 
all  directions. 

The  1/e  spatial  decorrelation  contour  for  these  two  propagation 
situations  is  shown  in  Figure  19.  For  reference,  two  radio  wavelengths 


FIGURE  19  PREDICTED  SPATIAL-DECORRELATION  PATTERN 


have  been  used:  one  la  the  observation  frequency  and  the  other  la  1000 
MHz.  For  the  observation  geometry  the  decorrelation  scale  la  large;  for 
the  higher  frequency  the  1/a  contour  Is  beyond  the  limits  of  the  figure. 

At  244  MHz,  the  contour  is,  of  course,  smaller  and  exhibits  the  pattern 
anisotropy  expected  for  the  azimuth  of  observation. 

The  field-aligned,  near-vertical- geometry  cases  are  more  per¬ 
tinent  to  DFA  applications.  Propagation  along  the  extended  axis  of  the 
irregularities  symmetrizes  the  spatial  correlation  and  reduces  the  decor¬ 
relation  scale  at  244  MHz  to  about  100  m.  /t  1000  MHz,  that  dimension 
Is  about  500  m.  For  our  purposes  here,  it  is  adequate  to  say  that  sig¬ 
nificant  distortion  of  the  desired  DFA  beam  should  be  expected  If  the 
aperture  dimensions  at  244  MHz  or  1000  MHz  approach  these  decorrelation 
scales.  Specific  calculation  of  sldelobc  degradation  will  be  discussed 
in  .Section  VI. 

We  can  extrapolate  from  the  contours  in  Figure  19  to  other 
Irregularity  anisotropies  and  turbulence  strengths.  For  example,  L-shell- 
aligned  propagation  through  sheetlike  irregularities  will  produce  a 
spatial  correlation  surface  narrow  along  one  axis,  but  extended  along 
the  other.  The  narrow  dimension  will  have  a  decorrelation  scale  similar 
to  that  for  field-aligned  propagation  and  rodlike  irregularities.  As 
for  changes  in  turbulent  strength,  the  correlation  contour  narrows  in 
all  directions  by  roughly  a  factor  of  two  from  that  in  Figure  9,  for  a 
10-dB  increase  in  Cg .  In  the  auroral  zone,  systematic  increases  of  this 
magnitude  seem  unlikely,  however. 

3.  Polar-Cap  Anisotropy  and  Drifts 

The  magnetic  polar  cap  is  the  region  of  open  magnetic  field 
lines  that,  driven  by  the  solar  wind,  extends  into  the  tall  of  the  mag¬ 
netosphere.  It  is,  accordingly,  a  region  in  which  the  structure  and 
dynamics  are  closely  tied  to  the  state  of  the  interplanetary  magnetic 
field  (IMF)  and  thus  to  the  sun.  This  solar  control  is  conspicuous  if 
we  consider  the  scintillation  changes  observed  over  the  recent  sunspot 
cycle.  A  measurement  campaign  in  1975  Judged  the  region  benign  in  terms 
of  communication  disruption,  but  in  1979  scintillation  levels  throughout 
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'llio  polar-cap  region  Ih  InLorcntliiA  not  only  because  of  the 
never  l.ty  of  hc  I  ntl  1 1  at  I  on ,  but  It  In,  very  likely,  a  nource  region  for 
the  large-scale  electron-density  ntructuron  obnerved  In  the  auroral  zone*. 
Furthermore,  the  dynamic*  of  the  region  provide  a  sensitive  Indicator  of 
the  sense  and  magnitude  of  the  Interplanetary  magnetic  field.  These  fac¬ 
tor*  are  essential  to  an  understanding  of  the  Irregularity  physics  and 
to  the  eventual  goal  of  global  l)FA  modeling. 

Until  the  observation*  that  are  discussed  here  wore  mado, 
virtually  no  measurements  of  Intermediate-scale  irregularity  anisotropy 
In  the  polar  cap  existed.  The  occasional  exceptions  are  from  the  early 
portions  of  Wideband  satellite  passes  In  the  midnight  sector  and  under 
disturbed  magnetic  conditions;  however,  the  1  ow-e levotl.on-anglo-observatlon 
geometry  for  these  Wideband  cases  Is  far  from  Ideal  for  the  detailed 
diagnosis  of  Irregularity  shape.  In  contrast  with  this  lack  of  anisotropy 
data,  the  plasma  drift  In  the  polar  cap  has  been  extensively  studied 
using  both  modeling  and  observation.  Experimental  concentration  has 
been  on  the  large-scale  electric-field  patterns  mapped  down  from  the 
magnetosphere.  These  have  been  Implied  from  multiple  magnetometer  data 
[ Fr I ls-Chr Is tensen  ct  al . ,  1972]  or  directly  measured  by  balloon  flights 
[Mozcr  and  Seri  In,  1969]  and  satellites  [Heel,  la  ct  a  l .  ,  1981]. 

The  polar-cap  plasma  flow  direction  is  dominated  by  the  north¬ 
ward  (B,,)  component  of  the  Interplanetary  field.  Under  normal  conditions, 

B y  Is  negative  anil  small,  this  results  In  a  convection  pattern  In  the 
vicinity  of  the  magnetic  pole  that  is  two-celled  with  flow  from  noon 
towards  midnight  and  sunward  return  through  the  dawn  and  dusk  auroral 
zones.  Generally,  the  two-flow  cells  are  asymmetric,  both  skewed  from 
the  noon-midnight  meridian,  and  offset  from  the  magnetic  pole.  This 
lack  of  symmetry  is  generated  by  spatial  or  temporal  changes  in  the  east¬ 
ward  component  of  the  IMF,  B^,  which  to  first  order  controls  the  dawn- 
dusk  flow  component. 

When  B.,  turns  northward,  a  variety  of  complex  clcctric-flcld 
distributions  can  occur,  lloppnor  fl972]  has  categorized  these  in  terms 


of  their  dawn-dusk  signatures.  These  range  from  patterns  explainable 
by  four-celled  polar-cap  convection  to  essentially  random  field  varia- 
tions  throughout  the  entire  region.  Burke  ct  al .  [1979]  have  studied  a 
subset  of  north  observations  and  have  postulated  a  magnetospher i< 
mechanism  for  their  cause. 

Satellite  in  situ  measurements  of  polar-cap  electric  fields 
such  as  those  referenced  above  have  been  made  using  sun- synchronous 
dawn-dusk  orbits.  These  provide  one-dimensional  snapshots  of  the  most 
critical  features  of  the  convection  dynamics,  but  cannot  provide  simul¬ 
taneous  information  about  the  field  elsewhere  in  the  polar  cap.  Tilt- 
power  of  the  ground-based  drift  measurements  using  the  polar  beacons  is 
that  the  drift  pattern  can  be  continuously  monitored,  albeit  at  a  single 
point.  The  complimentary  nature  of  the  two  methods  is  clear,  and  should 
be  taken  full  advantage  of  when  possible. 

The  data  we  discuss  here  are  from  an  extensive  spaced-rece i ver 
data  set  collected  from  Thule,  Greenland.  Since  1979,  short  spans  of 
multiple  antenna  observations  have  been  routinely  made  daily  for  several 
months  at  a  time;  our  emphasis  here  is  on  data  from  several  days  of  con¬ 
tinuous  observations  made  during  January  1982.  The  diurnal  drift  patterns 
derived  from  these  data,  using  signals  from  the  UHF  polar  beacons,  is 
currently  being  consolidated.  Here  we  present  anisotropy  and  drift  data 
for  one  consecutive  dawn,  noon,  and  dusk  sequence  during  which  the  mag¬ 
netic-field  conditions  were  nearly  constant. 

Dawn-sector  anisotropy  data  for  22  January  1982  is  shown  in 
Figure  20.  The  lower  scale  on  the  figure  is  universal  time  while  the 
upper  scale  indicates  corrected  magnetic  local  time.  Because  the  pene¬ 
tration  point  remains  nearly  constant  during  the  observations,  the  two 
time  scales  are  separated  by  a  fixed  amount.  The  invariant  latitude  at 
F-region  penetration  altitudes  is  near  85°,  well  into  the  central  polar 
cap. 

The  ground-diffraction-pattern  orientation  data  in  the  top 
portion  of  Figure  20  distinctly  indicates  rodlike  irregularities  as  the 
cause  of  the  scintillation.  In  fact,  the  model  predictions  cannot  be 
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FIGURE  20  LOCAL-DAWN-SECTOR  MEASURED  AND  MODEL  PATTERN 
ANISOTROPY 

distinguished  from  the  measured  data;  those  for  sheetlike  irregularities 
(not  shown)  differ  hy  nearly  30°  from  the  measured  orientation.  The 
irregularity  extension,  for  which  two  model  predictions  are  shown,  is 
quite  variable,  hut  generally  is  greater  than  that  seen  in  the  auroral 
zone  [ Livingston  et  al . ,  19821,  reaching  20:1:1  beyond  0330  UT.  The 
corresponding  dawn-sector  drifts  are  shown  in  Figure  21.  The  drift  mag¬ 
nitudes  are  high,  and  exceed  900  m/s  slightly  past  0500  UT.  The  magnetic 
north-south  component  systematically  increases  from  200  m/s  to  about 
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FIGURE  21  LOCAL-DAWN -SECTOR  MEASURED  IRREGULARITY 
DRIFTS 


500  m/s  southward,  while  the  more  rapidly  varying  east-west  component 
averages  about  400  m/s  in  an  eastward  direction.  For  this  invariant 
latitude,  the  southeastward  flow  is  consistent  with  an  overall  two-cell 
flow  pattern  expected  for  negative  conditions.  A  comparison  of 
Figures  20  and  21  suggests  the  same  convective  control  of  anisotropy 
observed  in  the  auroral  zone  [Livingston  et  al.,  1982].  Generally,  the 
axial  ratio  follows  the  variations  in  north-south  drift,  minimizing  when 
the  drifts  are  high.  The  association  between  anisotropy  and  the  east- 
west  velocity  component  is  not  clear.  This  point  will  be  discussed 
further  in  association  with  the  noon-  and  dusk-sector  data. 

The  local-noon-sector  data  correspond  to  nearly  the  same  invar¬ 
iant  latitude  of  85°  and  are  shown  in  Figures  22  and  23.  The  irregulari- 
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FIGURE  22  LOCAL-NOON-SECTOR  MEASURED  AND  MODEL  PATTERN  ANISOTROPY 


tics  are  still  rod  1 i kc ,  but  arc  extended  over  those  seen  in  the  dawn 
sector,  and  exceed  50:1:1  at  times.  An  in  situ  extension  of  20:1:1  is 
a  typical  value.  The  axial  ratio  varies  with  a  60-to-80-min  period  up 
through  1300  UT,  perhaps  associated  with  the  passage  of  F-region  sun- 
aligned  arcs  [Weber  and  Buchau,  1981],  The  drifts  in  Figure  23  are 
mostly  anti  sunward,  as  expected  in  the  noon  sector ,  but  also  include  a 
distinct,  but  rapidly  varying  westward  component.  This  presumably  arises 
from  an  asymmetry  in  the  two-celled  pattern  generated  by  a  nonzero  com¬ 
ponent  of  By.  The  magnitude  of  the  drifts  occasionally  exceeds  800  m/s, 
and  are,  therefore,  comparable  to  those  observed  in  the  dawn  sector. 
Unlike  the  dawn  sector,  however,  the  variations  in  axial  ratio  are  tied 
more  closely  to  changes  in  east-west  drift. 
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FIGURE  24  LOCAL-DUSK-SECTOR  MEASURED  AND  MODEL  PATTERN 
ANISOTROPY 


A  third  set  of  data  taken  near  local  dusk  on  the  same  day  is 
shown  in  Figures  24  and  25.  The  penetration  latitude  in  this  case  is  86 
slightly  northward  of  the  ottier  time-sector  data.  The  irregularities  ar 
still  rodlike,  systematically  increasing  from  about  20:1:1  elongation  in 
m id-af t ernoon  to  about  40:1:1  by  2100  UT.  This  increase  is  associated 
with  a  general  decrease  in  drift  speed,  although  there  is  no  clear  asso¬ 
ciation,  feature-for-feature  with  the  northward  or  eastward  components. 
The  northeastward  drift  in  this  case  also  fits  the  skewed  two-cell  flow 
pattern  suggested  by  the  dawn-  and  noon-sector  data. 
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FIGURE  25  LOCAl-DUSK-SECTOR  MEASURED  IRREGULARITY 
DRIFTS 


Detailed  comparisons  between  these  data  and  simultaneous  iono- 
sonde  and  optical  observations  are  underway.  Associations  between  partic 
ular  anisotropy  and  drift  features,  and  the  large-scale  polar-cap  struc¬ 
tures  can  then  be  established  and  the  process  of  identifying  instability 
processes  can  be  initiated.  From  these  measurements  alone,  however,  we 
can  make  some  preliminary  conclusions  appropriate  to  the  DFA  problem: 

(1)  The  irregularities  in  the  'V  600  m  spatial  wavelength 
regime  are  almost  exclusively  rodlike  and  range  in 
extension  from  about  1:1:1  to  10:1:1.  We  note  that 
on  the  day  following  our  examples  here,  deviations 
from  the  rodlike  structures  do  occur  correspond i ng 
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to  unusual  drift  flow  patterns.  For  either  irregu¬ 
larity  type,  scintillation  levels  for  field-aligned 
propagation  can  be  greatly  and  consistently  enhanced. 

For  the  polar  cap  this  will  he  at  a  near-vertical 
geometry. 

(2)  The  measured  drift  patterns  in  the  central  polar  cap 
fit  the  expected  two-celled  pattern  and  can  exceed 
800  m/s  in  the  dawn,  noon,  and  dusk  sectors. 

(3)  An  association  exists  between  drift  magnitude  and 
direction,  and  irregularity  anisotropy. 

Convection  dominance  of  anisotropy  has  been  clearly  demonstrated  in  the 
auroral  zone,  and  presumably  a  similar  process  is  ongoing  in  the  polar 
cap.  The  large-scale  dynamics  of  high-latitude  fields  are  relatively 
well  understood  and  are  effectively  modeled.  As  pointed  out  in  Section 
B,  if  the  anisotropy /dr i ft  associations  can  be  defined,  it  will  greatly 
simplify  the  global  prediction  of  irregularity  occurrence  and  DFA  per¬ 
formance  . 

4.  Polar-Cap  Turbulence  Strength 

Soon  after  the  initial  polar  beacon  observations  from  Thule  in 
late  1979,  it  was  clear  that  the  turbulence  levels  were  routinely  exceed¬ 
ing  those  simultaneously  observed  in  the  auroral  zone.  To  illustrate  that 
difference  here,  and  to  provide  direct  estimates  of  DFA  beam  distortion 
in  the  polar  cap,  turbulence  levels  for  the  noon  sector  on  Day  22  have 
been  calculated  for  use  with  the  anisotropy  and  drift  data  in  Figures 
22  and  23. 

The  measured  intensity  scintillation  index  and  calculated  C 

s 

levels  are  shown  in  Figure  26.  The  S.  index  and  C  trends  are  generally 

4  S 

well  correlated,  although  C^  exhibits  a  wider  dynamic  range  than  the 

index,  which  is  saturated  much  of  the  time.  The  peak  C  levels  here  are 

21  s 

near  10  ,  significantly  stronger  than  those  observed  in  the  auroral  zone, 

but  still  weaker  than  peak  values  seen  at  the  equator  [Livingston  et  al., 
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We  again  consider  two  propagation  situations,  and  the  correspond¬ 
ing  1/e  spatial  correlation  contours  are  shown  in  Figure  27.  The  first 
of  these  is  for  the  average  20:1:1  measured  axial  ratio  and  observation 
geometry.  The  correlation  surface  in  this  case  is  highly  anisotropic, 
with  a  minimum  l/e  width  of  about  200  m  at  244  MHz.  At  1000  MHz,  the 
pattern  would  be  of  similar  shape,  but  wider.  The  second  case  in 
Figure  27  is  for  field-aligned  propagation  and  more  extended  (50:1:1) 
irregularities.  At  244  MHz,  the  degradation  of  performance  would  be 
significant  even  for  small  DFA  apertures.  For  the  1 000  MHz  case  shown, 
the  decorrelation  scale  is  about  200  m,  much  shorter  than  that  in  the 
auroral  zone,  but  still  not  a  severe  constraint  on  moderate  size  DFA 
systems . 

The  scaling  of  the  contours  in  Figure  27  for  different  C 

s 

levels  is  similar  to  that  for  the  aurora]  data,  i.e.,  a  10-dB  increase 

in  C  will  roughly  halve  the  decorrelation  scale.  The  irregularities 
s 

sampled  for  this  data  set  are  produced  in  the  polar  cusp  and  have  drifted 


FIGURE  27  PREDICTEO-SPATIAL-DECORRELATION  PATTERN 
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into  the  central  polar  cap.  In  the  polar  cusp  itself,  turbulent  levels 

21 

will  undoubtedly  exceed  the  C  =  10  -dB  level  assumed  here,  and  may  be 

23  S 

comparable  to  the  >10  -dB  levels  observed  at  the  equator.  Decorre¬ 
lation  scales  less  than  100  m  could  result,  making  the  cusp  a  problem 
region  for  even  moderate-sized  apertures.  Clearly,  observations  such 
as  those  discussed  here  should  be  made  directly  in  the  polar-cusp  region. 
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VI  SUMMARY  AND  DISCUSSION 


In  this  report,  we  have  identified  the  principal  scintillation 
parameters  that  affect  the  performance  of  discrete  focused  arrays.  We 
then  reviewed  recently-analyzed  auroral-zone  and  polar-cap  scintillation 
data  to  determine  the  potential  for  DFA  performance  degradation  in  those 
regions.  We  showed  that  the  equivalent  isotropic  turbulence  strength 
Cg,  together  with  the  power-law  index  parameter,  sets  the  level  of  the 
perturbation.  The  irregularity  anisotropy  and  the  propagation  geometry 
dictate  the  final  degradation  of  the  DFA  performance. 

For  a  nonadaptive  system,  the  time  structure  of  the  scintillation 
affects  only  the  detectability  of  the  target.  In  an  adaptive  system, 
the  spatial  characteristics  are  varied  in  response  to  the  current  signal 
inputs.  In  that  case,  the  time  structure  of  the  scintillation  is  impor¬ 
tant,  but  we  have  not  addressed  adaptive  systems  in  this  report. 

In  Sections  V-B-2.  and  V-B-A.  we  show  that  moderate  to  severe  pro¬ 
pagation  disturbances  correspond  to  equivalent  turbulence  levels  over  a 
20  21 

100-km  layer  of  10  and  10  respectively.  Figures  19  and  27  show  the 
two-dimensional  spatial-coherence  structure.  To  determine  the  resultant 
degradation  of  DFA  performance,  we  have  used  Equation  (17)  for  a  linear 
DFA  aligned  along  the  axis  of  maximum  spatial  coherence  reduction.  For 
a  uniformly  weighted  array, 

N-l 

<|VR(mAK)|2>  =  N  ^  1  -  |^|  exp|-^C2<<)[n6£|P"1|exp  2m.^  (39) 

n=- (N-l ) 

where  A£AK  *  2*/N.  If  C?,  =  0,  Eq.  (39)  reduces  to  Eq.  (16). 

0<J> 

Figure  28  shows  the  pattern  degradation  of  a  16-element  array  with 
A£  »  10  m.  To  calibrate  the  spatial  wave-number  scale,  consider  that 
the  nulls  occur  at  multiples  of  x/NA£.  The  propagation  effects  broaden 
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Che  main  beam  and  reduce  the  depth  of  the  nulls.  Thi s  is  the  average 
degradation,  however,  and  individual  realizations  will  show  significant 
variations.  In  any  case,  severe  performance  degradation  would  he  expected 
for  the  moderate-to-severe  perturbation  levels. 

The  current  state  of  scintillation  modeling  and  the  available  data 
base  does  not  permit  the  actual  probability  of  occurrence  of  severe  DFA 
performance  degradation  to  be  specified.  We  do  know,  however,  that 
severely  disturbed  regions  are  confined  in  latitude  and  strongly  affected 
by  convection  and  a  conducting  E-layer.  The  rudimentary  elements  of  a 
global  structure  model  have  been  described  by  Kelley  et  al .  [1982]. 

The  planned  DNA  HILAT  satellite  experiment  should  provide  the  necessary 
data  base  for  implementing  such  a  model. 

We  also  note  that  a  single  power-law  model  does  not  accurately 
reflect  the  in  situ  irregularity  structure;  rather,  the  spectrum  has  a 
two-component  structure  with  shallowly  sloped  and  steeply  sloped  regiments 
separated  at  a  wave=number  corresponding  to  scale  sizes  between  500  m 
and  1  km  [ Rino  et  al . ,  1981;  Kelley  et  a  1  . ,  1982;  Basu  e t  al . ,  1982]. 
Accommodating  the  two-component  structure  does  not  affect  the  spatial 
coherence  measures  significantly,  however,  and  the  single  component 
model  can  be  used  effectively  for  deriving  DFA  system  performance  measures. 

We  suspect  that  an  understanding  of  the  detailed  spectral  character¬ 
istics  will  ultimately  improve  our  ability  to  predict  the  occurrence  and 
lifetime  of  structure.  Significant  improvements  in  our  predictive  capa¬ 
bility  must  await  such  refinements. 
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